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■ constants with the mass of the meson, as seen in the available experimental data. We use an 

Q-ii effective light-front QCD-inspired dynamical model regulated at short-distances to describe the 

valence component of the pseudoscalar mesons. The experimentally known values of the mass 



splittings, decay constants (from global lattice-QCD averages) and the pion charge form factor up 



^ . to 4 [GeV/c] are reasonably described by the model. 
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I. INTRODUCTION 



The hadronic structure viewed through effective hght-front theories inspired by Quan- 
tum Chromodynamics can shed hght in the investigation of the interaction 
between the hadron constituents and in the study of the transition from effective to funda- 
mental degrees of freedom that should be revealed at large momentum scales. The effective 
QCD-theory is expressed through a squared mass operator acting on the valence component 
of the hadron light-front wave function. The effective interaction embeds, in principle, all 
the complexity of QCD through the coupling of the valence state with higher Fock-states, 
reduced to the valence sector [5|. The effective squared mass operator acting depend on few 
physical parameters: the constituent quark masses, the effective quark-gluon coupling enter- 
ing in the Coulomb-like mteraction and the strength of the short range hyperfine interaction, 
fixed by the pion mass |6j . Hadronic observables are functions of these parameters, and one 
physically sensible is the dependence on the quark masses, which allows to get insight on 
the limit of heavy quarks. In particular, this limit is reflected in the weak decay constant of 
pseudoscalar heavy-light mesons, which in potential models were found to be oc 1/ ^ymQ |3] 
{rriQ is the heavy quark mass). Also, this scaling property has been shown in a light- front 

in our references.) 

The weak decay constant is closely related to the physics at small distances contained in 
the valence component of the pseudoscalar meson light-front wave function. It constitutes 
and important source of information on the short-range part of interaction between the 
quark and antiquark. Exjjerimental values of the weak decay constant are known for the 
pion, kaon, and D'^ The few parameter dependence of the effective theory can be 
translated into correlations between observables of a particular hadron or among different 
hadrons. Therefore, it is possible to indicate relevant relations between physical quantities 
that otherwise would have no simple reason to exhibit a close dependence, besides being 
properties of the same basic theory. For example, in the recent review of the application 
of Dyson-Schwinger equations to QCD, it was shown systematic correlations between differ- 
ent meson properties with mass scales, which was also useful to compare the pseudoscalar 
decay constants with results from Lattice QCD. 

The experimental values of the mass splitting between the ground states of pseudoscalars 
and vector mesons presents a systematic dependence on the corresponding pseudoscalar 
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mass, which is well described by the effective QCD-theory even without confinement p. 
The mass splitting is associated with the binding energy of the quark-antiquark pair in the 
meson, as that model lack a confining interaction. (It is worthwhile to note that the model 
account for the binding energy of the spin 1/2 ground state baryons containing two light 
and one heavy quark[ll|). 

In Ref. [12], its was pointed out that fps should scale with the sum of the constituent 



■ight-front QCD-inspired model it was 
13| |. In the light-front QCD-inspired 



quark masses, and more recently in the context of a 
found that fps scales with the vector meson mass 
model without a short-range regulator it was assumed the dominance of the asymptotic 
wave function Both frameworks assume that the log- type singularity in the matrix 

element of the axial current between the vacuum and the meson state is fixed by the pion 
decay constant /tt- In these approaches, the weak decay constant depends directly on the 
constituent masses and on the short distance component of the valence part of the light- 
front meson wave function parameterized by fj^. Although, the experimental results for 
light mesons up to D [9] suggests such increase, relativistic constituent quark models in 



the heavy-quark limit (see Refs. [A and numerical simulations with quenched lattice- 
QCD 14| indicate that > /^Jli^, which is still maintained with two flavor sea quarks 
|l4 llSj l and in the most recent global average of lattice results This behavior of the 
weak decay constants of the heavy-light pseudoscalars is also found in a Dyson-Schwinger 
formalism applied to QCD, where general arguments says that in the heavy quark limit 



fps should be cx 1/ y AL,|l7l| {Mps is the pseudoscalar mass). In order to study the mass 
dependence of the weak decay constant we can attempt to use a regulated form of the 
light-front constituent quark QCD-inspired model 3, y, |6|. In this case, the systematical 
investigation of the mass dependence of meson observables, can be easily performed as the 
masses of constituent quarks acts as model parameters, which can varied while the effective 
quark-gluon coupling entering in the Coulomb-like interaction is flavor independent. 

The short distance interaction between the constituent quarks in the squared mass op- 
erator equation of the effective light-front QCD-theory if regulated allows a finite result 
for the decay constant and electromagnetic form factor. In this case, the eigenstate of the 
effective mass operator, i.e., the valence component of the light-front wave function, would 
decrease faster than pj^ for large transverse momentum, which is enough to make finite the 
one-loop integration in the weak decay constant and form factor. One can get some infor- 
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mation on the short- distance behavior of the valence component of the pseudoscalar meson 
from the electromagnetic form factor, which is experimentally well known for the pion (see 
Ref. (l^), while for the kaon data exists below 0.15 (GeV/c)^ jiol. l2ol|. In particular, when 
the asymptotic wave function is assumed for the soft-pion, its radius and decay constant are 
related by = v^/(27r^) HQ- 

Our aim in this work, is to study systematically the mass dependence of the pseudoscalar 
weak decay constant, electromagnetic form factor and the mass splitting between the ground 
states of pseudoscalar and vector mesons, using a light-front QCD inspired model regulated 
at short distances. We choose the regulator in a separable form to simplify the formalism. 
The effective mass operator equation for the valence component of a constituent quark- 
antiquark bound system was derived in the effective one-gluon-exchange interaction approx- 



imation fll] and simplified in Refs. ^, |^ . The squared mass operator includes a Coulomb-like 
and a Dirac-delta hyperfine interactions acting on the spin singlet state responsible for the 
mass separation between pseudoscalar and the vector meson states. Here, we extend the 
model by introducing a regulator in a separable form in the singular part of the interaction. 
Then, the eigenvalue equation for the effective squared mass operator is written as: 

,.2,/ r^ ,.2,/ r^ f dx' dk'.O (x')9(l - x') 



( Amim2 a x /„^2\ 



\Jx{l — x)x'{\ — x') 
\g{Ml)g{M'Ai,{x'A) , (1) 



V 37r2 Q2 

mi and m2 are the constituent quark masses. The free squared mass operator in the meson 
rest frame is 

X 1 — X 

and Mg2 has primed momentum arguments. The form factor of the separable regulator 
function is qIMq). The projection of the light-front wave-function in the quark-antiquark 
Fock-state is given by tp. The mean four-momentum transfer is Q^. The strength of the 
Coulomb-like potential is proportional to a and the coupling constant of the regulated Dirac- 
delta hyperfine interaction is given bvA. Note that for 5'(Mq) = 1, the original unregulated 
form of the model presented in Refs. 0,0] is retrieved. (In Refs. and Q were used a local 
Yukawa potential for the regularization of the contact interaction, here we use a separable 
form for simplicity.) 



The dependence of the form factor in terms of Mq appears to be natural in a hght- 
front theory, in which the virtuahty of an intermediate state is measured by the value 
of the corresponding free squared mass. In the rest frame of the quark- ant iquark pair 
Mq = Pq" P+ and therefore proportional to the free value of Pq" - the minus component of 
the free momentum {P^ = Pq ± Pq ). 

Although, it is known a more developed form of the model which contains the explicit 
confinement we will be content in solving Eq. which is enough for our purpose of 
studying only the ground state. In practice from the solution of Eq. ((T)), the constituents 
quarks are bound and, in that sense, confined in the interior of the mesons ^. 

The present light-front model is a drastic approximation to a severe truncation of the 
Fock-space in the effective theory. In the initial truncation of QCD only one-gluon exchange 
was kept, which includes Fock states with up to qq plus one gluon, leaving out the complex 

spin-dependence and momentum-dependence in the 
hyperfine interaction are greatly simplified to get Eq.(^ (with g{M^)) = 1) and confinement 
is absent in the model. Therefore, the success of model should be understood as an useful 
guide in the investigation of mesonic properties which present a systematic behavior that 
depends only on few basic quantities, that are parameterized in the effective theory. 

This work is organized as follows. In sec. II, the QCD-inspired model is transformed to 
the instant form representation and the eigenvalue equation for the squared mass operator 
is solved. The valence component of the meson wave function is derived. In sec. Ill, we give 
the formulae for the electromagnetic form factor and weak decay constant, derived from 
an effective pseudoscalar Lagrangian used to construct the spin part of the pseudoscalar 
meson wave-function. In sec. IV, we present and discuss the results obtained with the 
regularized model for the mass splittings between the pseudoscalar and vector mesons, the 
weak decay constants and the pion and kaon form factors. Also, in sec. IV we summarize 
our conclusions. 



II. THE QCD-INSPIRED MODEL IN INSTANT FORM REPRESENTATION 

The effective mass operator equation for the lowest Light-Front Fock-state component of 
a bound system of a constituent quark and antiquark is rewritten in terms of the instant 
form momentum. Here we follow closely Ref. j6|. The general transformation from the 
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light-cone momentum to three- momentum was derived in Ref. The form of Eq. (^J) in 
the instant form momentum basis is particularly simple and convenient for the numerical 
solution, when the momentum carried by the effective gluon is approximated by a rotational 
invariant form. The momentum fraction is transformed to 

with unchanged. The individual energies are Ei = \Jm1 + k"^ (i=l,2) and k = \k\. The 
Jacobian of the transformation (x, k_[_) to k is: 

dxdk± = ,,f} dk , (4) 
mrA[k) 

where the dimensionless phase-space function is 

A(k) = ; (5) 

and the reduced mass is = mim2/(mi + m2). 

Using the momentum transformation defined above, the eigenvalue equation, (^, written 
in the instant form momentum basis is: 

' ' V^t' ^A(k)A(k')Q^ 7nr^A(k)A(k') ) 

where rrts = mi + m2, Mq = Ei + E2 and Mq has the primed momentum arguments. 
The square momentum transfer is approximated by the rotationally invariant form = 



— \k — fc'p. The phase-space factor is included in the factor 1/ y A{k)A{k') 
The valence component of the light-front wave function is 



ip{x, k±) 



^^'^ m. (7) 



^ x(l — X, 

The higher Fock-state components of the light-front wave function of the composite system 
can be expressed in terms of the lower ones, as shown by the method of the iterated re- 



solvents 14] (presented in greater detail in 



) and by a quasi-potential expansion on the 



light-front of the Bethe-Salpeter equation j25|. Therefore, it is possible to reconstruct re- 
cursively all the Fock-state components of the wave function from the valence component. 
In this way, the full complexity of a quantum field theory can in principle be described by a 
light-front effective Hamiltonian acting in the lowest Fock-state component of a composite 
system. 



A. Meson Valence Wave Function 



To easily manipulate and solve Eq. (jH)), it is convenient to work with the operator repre- 
sentation: 



(m^ + V + V^) \if) = M^\^) . 



The matrix elements of the Coulomb-like potential V are given by: 

{k\V\k') = ^""-^ " 



and the for the short-range regularized singular interaction one has: 



(8) 



(9) 



(10) 



Just for convenience we kept the the same superscript 5 in the short-range part of the 
interaction as in Ref. Q], although it is regulated here. We introduce a form factor defined 



by (^Ix) = di^o) / y ^(f^)^ which now includes the regulator. 

The eigenstate of the squared mass operator, (jH)), is trivially given by: 



(11) 



where G^{M. 



psJ 



M2 

ps 



Ml - V 



is the resolvent of the operator Mq + V . The charac- 



teristic equation for the eigenvalue of the squared mass operator is: 

A-i = ±(;^|g^(m2 = — Idkl dk' ii:-\ry(M^ M0\ EiMl' 

rrir ^ nir 



:i2) 



'm ^A{k') 

We have not yet defined A in Eq. p2|l . To do that, we first remind the characteristic 
equation of the renormalized theory with the singular hyperfine interaction ((^(Mq) = 1). 
The bare coupling constant is obtained from the value of the pion mass and substituted in the 
characteristic equation which gives the mass of the pseudodscalars. Then, the characteristic 
equation appears in a subtracted form, in which the divergence in the momentum integration 
is removed ^: 



1 
1 

rrir 



dk / dk' 



dk I dk 



m 
1 



{^G"" {Ml)\k') 



{k\G\Ml)\k') 



A{k') 
1 



{mu,m—) 



(mi,m2) 



(13) 



where is the mass of the hght constituent quark. Observe that, the physical information 
contained in the pion wave function at short distances is carried to any other quark-antiquark 
system in Eq. (fT^ by the operator 



0.{Ml) 



m 



(14) 



(m„ ,m-) 

which has its matrix element evaluated at the origin in (|T!^ . The hat indicates the operator 
quality. 

In the case of the present regulated model, we define for each meson a value of A assuming 
that the form-factor qIMq) selects the relevant momentum region of the interacting quarks, 
or the relevant region of virtuality of the quark-antiquark pair, within the particular meson. 
Thus, the matrix element of the operator C(M^) should be taken between states defined by 
the function glM^). Introducing the operator (9ps(Mp^) for a general pseudoscalar meson, 
which has expression analogous to Eq. (fT^ . one has: 

1 



psJ 



1 



=G^(M 



psJ 



m. 



(15) 



(mi,m2) 

Then, using the operators defined in Eqs. ()14p and (fT^ . it is reasonable to generalize 
Eq. ()13|) to the following form: 



ps 



{g\0^{Ml) - 0,s{Ml)\g) 



I ps 







(16) 



where {k\g)ps := gi^M^) with Mq = yk"^ + mf + yk'^ + The strength of the short-range 
interaction for each pseudoscalar meson is determined by the pion mass and the regulator 
form-factor, according to: 



A;;= ps{gmM'^)\g) 



ps • 



(17) 



In the three-momentum basis Eq. (fT^ reads: 

/ 

dk [ dk' g{Ml) 



1 1 



{k\G\Ml)\k') 



rrir- 



4W) 



rrir 



A{k) 



IG^'iMDlk') 



1 



A{k') 



(18) 



(mi ,m2)/ 



where Mq and Mq are computed for the quarks with masses mi and 1712. 



In our calculation procedure, the resolvent is numerically obtained from 



C'iMl) = Go{Ml)+G,{Ml)T''{Ml)GM 



psl 1 



where the T-matrix is the solution of the Lippman-Schwinger equation: 



T^(M^2J = V + VGo{Ml)T^{Ml) 



psJ 



psJ 



(19) 



(20) 



where the free resolvent is Gq{M 



psJ 



ps 



The detailed expressions can be found 



in Ref. J 

The valence component of the light-front wave function of the meson is the solution of 
Eq. (dl) given by Eq. (fTT|) . which we write explicitly as, using Eq. (fT^ : 



1 



G 



ps 



X{1 - X) Ks - ^0 



g{Ml) + j dk'^ ^^{k\T''{Ml\k')gM') 



(21) 



where the overall normalization factor of the qq Fock-component of the meson wave-function 
is Gps- The three-momemtum is expressed in terms of the light-cone momentum with the 
transformation (jH)). The first term in Eq. fIT\i dominates for large momentum transfers if 
5f(Mo) = 1 (corresponding to the asymptotic form), differently from this situation when 
^^(Mq) 7^ 1 the two terms can compete even in the asymptotic region. 



B. Constituent Quark Masses and Mass Splittings 

Within the present model, the low-lying vector mesons are weakly bound systems of 
constituent quarks while the pseudo-scalars are strongly bound. Therefore in this model, 
the masses of the constituent quarks are obtained directly from the vector meson masses, 
as bJ: 



ruu 


2 ^ 


384 MeV , 




rris 


= Mk*- 


^Mp = 508 MeV , 




rric 


= Md^ - 


-M„ = 1623 MeV , 
2 ^ 




TJlb 


= Mb* - 


^Mp = 4941 MeV , 


(22) 



where the masses of the vector mesons are 768 MeV, 892 MeV, 2007 MeV and 5325 MeV 
for the p, K*, D* and B*, respectively Q]. The constituent masses for the up and down 
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quarks are considered equal (we disregarded the small few MeV difference in the current up 
and down masses j^). Using the value of the light-constituent quark mass of 384 MeV, and 
assuming that the effect of chiral symmetry breaking is about the same for each flavors one 



gets an estimate of the current quark mass as itIq'^"^ = mq — rriu 



The current quark 



masses are estimated as m™^^ = 124 MeV, m^"^^ = 1239 MeV and mf"^'' = 4557 MeV 

In our model, the binding energy of the constituent quarks in the pseudoscalar mesons, 
is interpreted as the ^Sq - ^Si meson mass splitting, and thus a quantity directly related to 
data. The binding energy is simply Bps = M^ — Mpg defined to be positive. The experimental 
values for the ground state quantities show evidence for a strong correlation of B^s and Mps 



qualitatively reproduced by the renormalized model with singular interaction 



We will 



see in sec. IV, that Eq. fll6|) also provides a reasonable description of the mass splitting. 



III. ELECTROMAGNETIC FORM FACTOR AND WEAK DECAY CONSTANT 



To obtain the e 
gestion of Refs. [22, 



ectromagnetic form pseudoscalar decay constants, we follow the sug- 



26|. To construct such observables, one describe the coupling of the 
pseudoscalar meson field to the quark field, by an effective Lagrangian density with a pseudo- 
scalar coupling between the quark {qi{x) and q2{x)) and meson {^ps{x)) fields: 

'^effix) = -iTps^ps{x) qi{x)j^q2{x) + h.c. , (23) 

where Tps is a constant vertex. After the integration in the minus momentum component of 
the momentum integration of the one-loop amplitudes that define the electromagnetic form 
factor and weak decay constant, the asymptotic form of the wave function is substituted by 
the valence component of the model wave function. The integration in the minus momentum 



component eliminates the relative time between the quarks in the intermediate states 



Form Factor of Pseudoscalar Mesons 



The pseudoscalar meson electromagnetic form-factor is obtained from the impulse ap- 
proximation of the plus component of the current = + j^) in the Breit-frame with 
momentum transfer g"*" = and = —(f satisfying the Drell-Yan condition. The general 
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structure of the qq bound state forming the meson comes from the pseudoscalar couphng 
(j2Sl)- We use such spin structure in the computation of the photo-absorption amphtude in 
the impulse approximation (represented by a Feynman triangle diagram), which is written 
as: 



X Y 



(2^ 

;^ + nil 



tr 



It + 1712 



k"^ — 7712 + (^ ~ ~ "^1 + 



1^' + TUi 



-7 



+ [l--2], 



(24) 



{k — pY — 7n\ + ie 

where Fps{q^) is the electromagnetic form-factor and Cj is the quark charge. The meson 
momentum in the initial and final states are defined by p° = p'^ and p' ^ = —p± = 
Nc = 3 is the number of colors. 

The choice of the plus-component of the current is adequate in the case of the pseu- 
doscalars, because after the integration over k~ = k^ — k^ the suppression of the pair dia- 
gram is maximal for this component in the frame wher e = and just the wave- function 



26 



271 128[. In our model only the valence 



components contribute to the form-factor j2^ 
component is considered. Although, we compute the integration in the minus momentum 
component assuming a constant vertex, one can identify in the expression how the valence 
component of the wave-function correspondent to the non-constant vertex of Eq. (PT|) should 
be introduced. As the details of this derivation is by now standard, we present directly the 
final result: 



ei 

k~^(mi 



1 dx 



d^k, 



(27r)3 7o l-x 

m2)^ - k^<f^) 



2mim2 - 27n\ + k^^^p'^ + 
tfjpsix, K±)tlJps{x, ifi) + [1^2] 



(25) 



where the momentum fraction is x = k^ /p'^ and k^^^ = {k\ + 7nf)/k^ . The quark transverse 
momentum in the meson rest frame is given by: 



K± = k± + X 



(26) 



and K'j_ = Kj_ — xq±. The expression for form factor gives the standard Drell-Yan for- 
mula once the bound-state wave-function of the constant vertex model (asymptotic form) is 
recognized 

Gr 



'ps 



^x{l -x{7nl- M2) 



(27) 
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which is the first term in Eq. ()21|) for 5'(Mq) = 1. The second term in ()21|) comes from the 
Coulomb-hke interaction. The other factors in Eq.()25p compose the Melosh rotations of the 
individual spin wave function of the quarks. 

The size of the meson is closely related to the square-root mean square charge radius 
which is calculated as 



4k) 



In sec. IV we adjust the regularization parameter (see Eq.()34|) by fitting the pion charge 
radius, \J (r^), which has the experimental value of 0.67 ± 0.02 fm 2^. The charge radius 
from Eq.(j2ni) in the soft-pion limit (m^ = 0) usin g th e asymptotic wave-function (jTfj) . with 



= V2 mu(d)/ f-K from the Goldberger-Treiman 



the well known expression Mrl) = V3/(27r/^) from Ref. j2jj 
factor pHjl for = 0, reduces to expression for as given in 
model, for a = 0, g{Ml) = 1 and = recovers the soft-pion result, i.e, J (r^) 



relation at the quark level, results in 
n this case, also the form 
We observe that our 
0.58 fm. 



B. Weak Decay Constant 

The leptonic weak decay constant of the pseudoscalar meson (/ps) is a physical quantity 
that depends directly on the probability to find the quark- ant iquark Fock-state component 
in the meson wave-function Also, fps depends on the short-range physics carried by the 
wave function when the quark and antiquark are close. 

The meson weak decay constant is calculated from the matrix element of the axial current 
A^(0) between the vacuum state |0), and the meson state \p)ps with four momentum p 0]: 

(0 I A^{Q) I = iV2UsT^ , (28) 

where A^{x) = tq{x)'j^'~f^q{x). 

The matrix element of the plus component of the axial current is derived from the pseu- 
doscalar Lagrangian, (PHj) . and it is expressed by a one-loop diagram, which is given by: 



f -NT f Tr[7V(^-|^ + m2)7^(|^ + mi)] 

IVMVlpshs l^clpsj (2vr)4((fc-p)2-mi+^£)(P-m2 + Z£)' 

the plus component is used to eliminate the instantaneous terms of the Dirac propagator. 

By integration of Eq. ()29|) over k~ , one obtains the expression of fps suitable for the 
introduction of the meson light front wave-function. So, performing the Dirac algebra and 
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separating the poles in the k -plane and integrating, one gets: 



Us = - ^N^l ^^Y^ ((1 - x)m, + xm,) J d k^J^^^ ^ (30) 

where quark 1 has momentum fraction x. This expression is written in the meson rest-frame 
and we have used the momentum fraction x = /p'^. The free square mass is defined in 
Eq. Note that this expression has a log- type divergence in the transverse momentum 
integration which was discussed in Ref. and parameterized in terms of /,r- 

One can write Eq. (jHn|) in terms of the valence component of the pseudoscalar meson wave 
function from Eq. ()2H1 . as: 

Us = I I ((1 - x)mi + xm2) j d^k^tp^x, k±) , (31) 



with X being the momentum fraction of quark 1. If one chooses g{M^) decaying as Mq^ for 
any r/ > is enough to make fps finite. 

In the particular case of 5'(Mq) = 1 the meson wave function, Eq. ()21|l . for large trans- 
verse momentum behaves as the asymptotic wave-function, which decreases slowly as pj^ 
implying in logarithmic divergences in the transverse momentum integrations of the weak 



decay constant and form factor. In Refs. 



12l | and jl^, the log-type divergent factors in 



the pseudoscalar decay constants were parameterized in terms of /^r and the sum of the 
constituent quark masses, which in the QCD-inspired model could be identified with the 
ground state vector meson mass. Therefore, one has: 

fps = const. / dx ((1 — x)m2 + xmi) , (32) 
Jo 

and oc + Q. In ou. n^odel ,n, + . the vecto. o,e=on o,a., and then as 

suggested by Ref. the decay constant scales as: 

(33) 



fps _ 



f. M, 

which approximates the existing data up to the kaon and D mesons but is not supported 
by relativistic constituent quark potential models and Lattice-QCD calculations as we have 
discussed in the introduction. The use of the separable regulator in the model brings this 
consistency. We will return to this discussion in the next section. 
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IV. DISCUSSION OF THE NUMERICAL RESULTS AND CONCLUSION 



The present QCD-inspired model for the effective squared mass operator acting only 
on the valence component of the meson wave function, Eq. ^ has the canonical number of 
parameters (the quark masses and a) plus two, when we choose the regulator form factor as 

The form factor (a) for equal mass quarks is the familiar Yukawa form in coordinate space. 
The other choice just contains the first two terms of a Taylor expansion of g^""^ for large Mq . 
In the limit of Z?^"-* — > oo the model with form factor g^^^ reduces to the renormalized model 
of Ref. 0], while by construction the form factor is qualitatively different as it does not 
allow to match that model. Nonetheless, for finite /5's as we will see from our numerical 
calculations both form factors produce practically the same results. 

The parameters P and the strength of the separable interaction \ps, Eq. (IT7|) . are ad- 
justed to reproduce the experimental pion charge radius 0.67 ± 0.02 fm [29] and mass, 
M7r= 140 MeV, for each form factor g^^^ and g^^^ independently. We use the light quark 
constituent mass of 384 MeV (see Eq. ()22|l ) and fj^ results 110 MeV for both separable 
interactions (see Table I) . The somewhat higher value for f^^ compared to the experimental 
value of 92.4±.07±0.25 (sl, is a common shortcomin g of light-front models of the pion when 
the valence wave function is normalized to one 



26, 



28|. The valence component appears 



to have a probability of about 70% (see also |26l. bsl]) which brings the model results for 
/,r to 92 MeV. In the case of a = 0.5, the resulting parameters for the form factors are 
(ji^) =-(634.5 MeV)2 and (3'-''^ =(1171 MeV)^. 

In figure 1, we show the results for the ^^i-^^o meson mass splitting, the binding energy 
Bps, as a function of the pseudoscalar meson mass. We choose a = 0.5 and the form 
factor regulator g^""^ (The differences between the masses obtained with the two form of 
regulators are less than 1 MeV). In the figure, we see the dependence of the mass splittings 
of qQ mesons with the pseudoscalar mass, obtained by the the variation of mg, while rriq is 
fixed at the values of 384 MeV (solid line), 508 MeV (dashed line) and 1623 MeV (dotted 
line). In this way, we simulate the families of mesons with an up or down, a strange and 
charm quarks and a distinct one which has mass rriQ. First, we compare the results of the 
present regularized model with the previous results of the renormalized model found for 
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TUq = 384 MeV and a = 0.5. In this case, the regularization increases Bps as seen in the 
figure. The regularization procedure naturally softens at short-distances the attractive part 
of interaction, which should be compensated by an effective increase of the strength of the 
separable interaction to keep the pion still strongly bound at its physical mass. The increase 
of the strength is reflected in the increase of binding, as seen in the flgure. Still the trend of 

he experimental values of the mass splitting for p — vr, K* — K^, D*^ — and B* — B^ 

^ is found. 

The results for the mass splitting for mesons containing at least one strange meson (dashed 
line in flgure 1) exhibit the same qualitative behavior found for mesons with an up or down 
quarks, i.e., the mass splitting decreases with the rise of the mass of the heavy quark. 
This should be the case since the masses of the constituents up-down and strange are very 
much similar, with an expected increase in the mass splitting when the up-down quark 
is exchanged with one strange quark which is heavier. By rising the mass of one of the 
constituents for mesons with charm (dotted line of flgure 1) the splitting increases, because 
the quarks become spatially closer and the binding is expected to rise as in nonrelativistic 
potential models. Also, as expected, the saturation of the binding energy appears for large 
masses. 

In flgure 2, we show the weak decay constant as a function of the intensity parameter a of 
the Coulomb-like interaction for different mesons. The calculations are performed with the 
regulator form-factor g^^^M^) = (6 + Mq)^^ with the parameter b adjusted for each given 
a between 0.1 and 0.5 in order to reproduce /t^ = 110 MeV. The kaon weak decay constant 
varies less than one MeV in this interval keeping the value 126 MeV (see Table I). We show 
in the flgure only results for (solid line), (dashed fine), B^ (solid line with dots) and 
B^ (dashed line with dots). The decay constants rise with a, as the qQ systems become 
more bound and compact due to increase in the Coulomb-like interaction. The effect is 
particularly dramatic for the heavier mesons B^ and B^ as could be anticipated thinking 
within a nonrelativistic potential model, where the probability to found the quarks at the 
origin should increase when the attractive force is strengthened. 

The pseudoscalar meson weak decay constant as a function of the vector meson mass 
is shown in flgure 3 as suggested by Eq. The calculations were performed with the 

regulator form-factor (a) and a = 0.5. In the flgure, we see dependence of fps for qQ mesons 
with the vector meson mass, obtained through the variation of rriQ, while niq is flxed at 
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the values of 384 MeV (solid line), 508 MeV (dashed line) and 1623 MeV (dotted line). 
The naive model of Eq. (jHHj) . which presents a linear increase of fps with M^, although it 
gives some qualitative insight into the data fails to describe the saturation and decrease of 
the results of the regulated model, which as expected has a fps decreasing with the mass 
of the meson. The data for is indeed below the linear curve and consistent with the 
dashed curve calculated with the regulated model for sQ pseudoscalars. There are several 
experimental values for fo+ obtained by different collaborations as quoted in Table I. In 
figures 3 and 4 we just indicate the experimental result from j3l| . 

The values of fps for mesons with one charm quark (dotted line in figure 3) increase with 
My, as the system becomes more compact up to the point that fps saturates for » 
mc= 1623 MeV (the probability density at the origin does not change anymore) while the 



expected 1/ yMg dependence dominates for large values of M^. 

In figure 4, the weak decay constant as a function of the pseudoscalar meson mass obtained 
in our regulated model with form factor (a) and a = 0.5 is compared to the recent global 
average of lattice-QCD results The short-dashed line gives a least-square fit to the 

experimental values of f-j^ and fx together with the lattice estimates for and Il4l 
given by /^^ = (0.0065 + 0.014Mp,)/(l + 0.055Mp, + 0.15M2J GeV^ as given in Ref. |lo|, 



where the 1 / y Mp^ behavior for large masses is built in. The results for uQ pseudoscalars are 
qualitatively agreement with that fit. Our calculations for uQ and sQ pseudoscalar mesons 
are in a good consistency with the global lattice averages of the weak decay constants, as 
seen by comparing the solid line with the full circles for uQ mesons and the dashed line with 
the full stars for sQ mesons. 

To close our study of the present regulated model in figures 5, 6 and 7 we show results 
for the pion and kaon electromagnetic form factors using a = 0.5 with the model regulated 
with form factor (a). The pion mean square radius is reasonable fitted and as well as the 
form factor up to about 4 [GeV/c]^ as shown in figure 5. The experimental values for kaon 
form factor |2^ present large errors and do not allow a definite conclusion as seen in 
figure 6. For completeness, we present the kaon form factor calculation up to 10 [GeV/c]^. 
We also compare with the calculations with the form factor (b), and we do not observe a 
strong model dependence below 4 [GeV/c]^. 

In Table I, we present the results for the pseudoscalar weak decay constants fps for tt, 
K, D^, D^, B^, Bg and B^ compared to global estimates of lattice-QCD results and ex- 
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perimental data. The consistence with lattice results indicates that the regularized model 
is able to parameterize the QCD-physics at short-distances in the ground state of the pseu- 
doscalar mesons quite reasonably. The pseudoscalar masses are underestimated for the heavy 
mesons, as seen already in figure 1, although the saturation behavior of the mass splitting 
that the data indicates is verified by the calculation. This problem can be overcome by the 



introduction of confinement in the model 



33|. 



In summary, we have shown that the suggested separable form to regulate the singular 
interaction in the square mass operator provides a reasonable description of the mass splitting 
between ^Si and ^5*0 meson ground states, the weak decay constants as found in a recent 
global average of lattice results [l^ and the pion form factor up to 4 [GeV/c]^. The main 
point here is that the model can describe the mass dependence of the weak decay constant, 
revealing that the physics in this observable is dominated by the mass of the meson itself, 
through the quark masses and binding. The effective squared mass operator acting on the 
valence component of the light-front meson wave function is again tested and proved to 
reasonably parameterize the dynamics of the constituents at short distances. The present 
version of the model does not have explicit confining interaction, therefore it is not able to 
account for the spectra. A more sophisticated version of the model that includes confinement 



was shown to describe the meson spectrum 



24| and the pion form-factor in the space and 



time-like regions [3^, can also be used in the future in a regularized form to allow the 
calculation of the pseudoscalar decay constants. 
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TABLES AND FIGURES 



QQ 


Aa) 
Jps 


Ah) 

Jps 


£ Latt n i^l 

Jps \m 


J ps 


Mps 


Mp/ [9J 




110 


110 




92.4 ± .07 ±0.25 [9] 


140 


140 


K+{us) 


126 


121 




113.0 ± 1.0±0.31[9j 


490 


494 


D^{cd) 


164 


159 


166±8±13li^(xlog) 


212il^Ili[9] 
202 ± 41 ± 17 [31] 
2621^^ ± 18 


1861 


1869 


Dtics) 


184 


178 


187±8±15 


188 ± 23 [2] 


1961 


1969 


B+{ub) 


118 


117 


135±16li0(xlog) 




5242 


5279 




154 


154 


156±18 




5342 


5370 


Bticb) 


375 


375 






6257 


6400±400 



TABLE I: Results for the pseudoscalar weak decay constants fpg compared with others calculations 
and experimental data. The calculations for fpg for a = 0.5 with regulator form-factors ^("•'(Mo) = 
(-(634.5)2 + M^y^ and g^^\Mo) = 1/M^ + (1171/M^)^ are given in the second and third columns, 
respectively. In the fourth column the global estimates from lattice-QCD results 1161 1 are shown. 
The experimental values for the decay constants 1311 . I32l | are shown in the fifth column. In 
the last two columns it is given the masses of the pseudoscalars obtained with model (a) and the 
corresponding experimental values ^. The masses obtained with model (b) (not shown in the 
table) and model (a) differ by less than 1 MeV. The decay constants and masses are all given in 
MeV. (The experimental errors in the masses are small excepting B'^ which has a large error.) 



20 



1000 I I 




^'^O 1000 2000 3000 4000 5000 

Mp, [MeV] 

FIG. 1: ^S'l-^S'o meson mass splitting (Bps) as a function of the pseudoscalar meson mass. 

The calculations are performed with the regulator form-factor g^^'^Mo) = (P^""^ + Mq)^^ (/3*^'*^=- 

(634.5 MeV)^ and a = 0.5) for qQ mesons and a varying mass mg with ruq fixed at 384 MeV 

(solid line), 508 MeV (dashed line) and 1623 MeV (dotted line). The results of the renormalized 

model from Ref. with a = 0.5 and fixed mg=384 MeV are shown by the short-dashed line. The 

'"I 

experimental values of the mass splitting for p — n, K* — K^, D*^ — and B* — B^ |3] are given 
by the full-circles. 
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FIG. 2: Weak decay constant as a function of the intensity parameter a of the Coulomb-hke 
interaction for different mesons. The calculations are performed with the regulator form-factor 
5(")(Mo) = + M§)-^ with the parameter fitted from f^, = 110 MeV for each given a. 
Results for (solid line), Df (dashed line), S+ (solid line with dots) and (dashed line with 
dots) . 
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FIG. 3: Pseudoscalar meson weak decay constant as a function of the vector meson mass. The 
calculations are performed with the regulator form-factor g(")(Mo) = {P^""^ + M^)~^ 
(634.5 MeV)^ and a = 0.5) for qQ mesons and a varying mass niQ with ruq fixed at 384 MeV 
(solid line), 508 MeV (dashed line) and 1623 MeV (dotted line). The short-dashed line give fps 



obtained from Eq. (|33|) . The experimental values are given by the full squares 
order of increasing values); the cross j31,] (/£>+) and the full triangle {fjj+)- 



(f-^ and fx in 
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FIG. 4: Pseudoscalar meson weak decay constant as a function of the pseudoscalar meson mass. 
The calculations are performed with the regulator form-factor q^^'^Mq) = (P^""^ + Mq)~^ (/3*^")=- 
(634.5 MeV)^ and a = 0.5) for qQ mesons and a varying mass mg with mq fixed at 384 MeV (solid 
line), 508 MeV (dashed line) and 1623 MeV (dotted line). The global estimates of lattice-QCD 
results |l^ for and /^+ are given by the full circles and for and /^o by the full stars. The 
short-dashed line gives a least-square fit to the experimental values of /tt and fx together with the 
lattice estimates for and B'^ 141 as performed in Ref. llOl l. The experimental values are given 

Q n 

by the full squares [9|] (/tt and fx in order of increasing values); the cross |31|] (/£>+) and the full 
triangle (/^+). 
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FIG. 5: Pion electromagnetic form factor. The results of the calculations performed with a=0.5 
considering the regulator form-factors c/^^H^o) = (-(634.5)^ + M|f)""^ and q'-^^Mq) = + 
(1171 /Mq)^ are given by the solid and dashed lines, respectively. The experimental data are from 
Ref. 11^. 
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FIG. 6: Kaon electromagnetic form factor. The results of the calculations performed with a=0.5 
considering the regulator form-factors c/^^H^^o) = (-(634.5)^ + M^)~^ and g^^\Mo) = + 
(I171/Mg)^ are given by the solid and dashed lines, resp ectively. The experimental data from 
Ref. 1^ are shown by the full triangles and from Ref. |20| by the full circles. 



26 



0.01 

2 4 6 8 10 

-q^[(GeV/cf] 

FIG. 7: Kaon electromagnetic form factor up to 10 (GeV/c)^. The results of the calculations 
performed with a=0.5 considering the regulator form-factors g^°'\MQ) = (—(634.5)^ + Mq) ^ and 
^W(Mo) = 1/Mq + (1171/Mq)^ are given by the solid and dashed lines, respectively. 
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